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INPUT ADMITTANCE OF A RECTANGULAR 
WAVEGUIDE -FED APERTURE ANTENNA RADIATING INTO AN 
INHOMOGENEOUS LOSSY DIELECTRIC SLAB 

By Calvin T. Swift 
Langley Research Center 

SUMMARY 

The input admittance of a rectangular waveguide opening onto a perfectly con- 
ducting ground plane and covered with an inhomogeneous dielectric slab is formulated as 
a boundary-value problem. As a specific application of the solution, a plasma slab is 
considered. In the calculations, the plasma is assumed to be homogeneous except in the 
vicinity of the ground plane, where a local boundary layer of electron density occurs. 
Two ranges of collision frequencies are selected, and calculations are performed for the 
input admittance as a function of peak electron density and boundary- layer thickness. 

The numerical results show that a small boundary- layer reduction of electron density 
near the ground plane can significantly influence the input admittance of the antenna and 
causes a substantial lowering of the voltage reflection coefficient when the plasm a is 
overd ense . 


INTRODUCTION 

The input admittance of a homogeneous plasma-covered rectangular waveguide-fed 
aperture antenna was originally investigated by Galejs (ref. 1) and Villeneuve (ref. 2) 
and then by Compton (ref. 3). Although calculations have been available for a number of 
years, only recently have attempts been made to confirm experimentally these theoreti- 
cal results. Using a Langmuir probe to infer electron density, Taylor (ref. 4) placed a 
rectangular waveguide- fed aperture antenna into a shock tube and measured the voltage 
reflection coefficient as a function of electron density. Best preliminary data for over- 
dense plasmas resulted in a measured reflection coefficient which was significantly 
lower than that predicted by homogeneous plasma calculations. Taylor concluded that 
the inconsistency was due to a local boundary- layer reduction of electron density at the 
ground plane. This conjecture was supported by the results of Meltz, Freyheit, and 
Lustig (ref. 5). They had previously observed a discrepancy between theory and experi- 
ment and attributed it to the presence of a thin glass wall between the cylindrical gap 
antenna and the plasma. 



These experimental observations suggested that some provision must be made in 
the th eory to account for variations of the dielectric constant near the antenna. This is 
done herein by considering the plasma as an inhomogeneous dielectric slab. The deriva- 
tion of the admittance expression closely follows the method given in reference 3. The 
fields are expressed in terms of their Fourier transforms, and the boundary-value 
problem is solved by assuming that the aperture field is the dominant waveguide mode. 
Parseval’s theorem is then used to obtain an admittance expression, consisting of a 
double integration over mode space, which must be numerically evaluated. The present 
formulation differs from that given in reference 3 in two respects. First, the fields are 
developed from single- component electric and magnetic vector potentials (instead of a 
two- component magnetic potential) and, second, the Helmholtz equation is numerically 
solved by utilizing the plane wave program of Swift and Evans (ref. 6). The numerical 
solution is required since the dielectric constant of the slab is assumed to be a compli- 
cated continuous function of one of the coordinates. 

Calculations are performed for the input admittance as a function of peak electron 
density and boundary- layer thickness for moderate and low values of the collision fre- 
quency. It is concluded that a thin boundary layer over the antenna can have a profound 
effect on the admittance and, indeed, reduces the reflection coefficient of the rectangular 
aperture covered with an overdense plasma. 

The special case of a homogeneous slab over the aperture is considered in the 
appendix. The excitation of surface waves in a nonlossy slab is discussed, and the 
admittance integrals are tabulated into forms suitable for numerical evaluation. 

SYMBOLS 


A electric vector potential 

A* magnetic vector potential 

a short dimension of rectangular waveguide 

b long dimension of rectangular waveguide 

bj n input susceptance ^normalized with respect to Yq^ 

c speed of light in free space 

E electric field intensity 
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frequency 


f(0,z),g(0,z) normalized Fourier transforms of vector potentials 

g in input conductance ^normalized with respect to Yq-^ 

g g surface-wave conductance ^normalized with respect to Y 0l) 

H magnetic field intensity 

i = fi 


c 0 


N 

N e 

N e ,o 

P 

V 0 

x.y,z 

Y in 

Y 0 

Yoi 


7 te 


wave number in free space, 


k x ,ky,k z Cartesian components of wave number 


complex index of refraction 

electron density 

peak value of electron density 

power 

aperture voltage 
Cartesian coordinates 
input admittance 


characteristic admittance of free space, 

characteristic admittance of the TEq^ mode of the rectangular waveguide, 
Yn 


1 - fa) 

kgbj 


I 

transverse electric component of admittance ^normalized with respect to 


C 01 
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y TM 

yin 

z o 

°P 

r 

€ 

X 0 

MO 

v 

<$> 

<t>* 

CO 

W P 

w p,o 


transverse magnetic component of admittance ^normalized with respect 

to *0l) 


input admittance (normalized with respect to 



thickness of dielectric slab 


attenuation constant of the plasma, 



polar components of k x /kg and ky/kg 
surface wave pole 
voltage reflection coefficient 
permittivity 

wavelength in free space, c/f 
permeability of free space 
collision frequency 
electric scalar potential 
magnetic scalar potential 
frequency of propagation 
plasma frequency, 2v x 8.97 x 10 
peak plasma frequency 
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Subscripts: 


I in dielectric slab 

0 in free space 

1 imaginary part 

r real part 

x,y,z vector component in direction indicated by subscript 

Superscripts: 

I vector component in dielectric slab 

II vector component in free space 

TE transverse electric 

TM transverse magnetic 

A double bar over a symbol indicates the double Fourier transform. A prime 
denotes the derivative with respect to the z -coordinate. An asterisk denotes complex 
conjugate except when used with magnetic potentials. 

THEORY 


The geometry of the problem is shown in figure 1. A rectangular waveguide, 
excited in the dominant TEqj mode, opens onto a perfectly conducting ground plane 
which is covered with an inhomogeneous dielectric slab. The dielectric constant of the 
slab is assumed to vary continuously in the z-direction for 0 I z S z 0 (region I), and 
air is assumed to fill the remaining half space for z g z Q (region n). 

It can be shown that the field components in both regions may be constructed from 
a set of vector potentials A^’ 11 and (a^/g) ’ , such that 
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where a time harmonic variation of the form e'3 wt is implicit and the vector potentials 
in region I satisfy the partial differential equations 


2 a I 1 d A I , „ 2 e I A I 


* z ej dz dz z 0 e 0 "z 


^ A z - -r^ -r- A, + k n * A„ = 0 

i*\I 


? 2( A *zY + k 2 6 I (z) f A zy_ 

+k o -ej-v— y " 




( 2 ) 


The vector potentials A z and A* are those defined by Stratton (ref . 7) . It should be 
noted that the scalar potentials cp and <p* are eliminated by the Lorentz conditions 
8A Z * 8 /A*\ 

jco^ 0 e <p = - 7 ^- and - — I— - J. The structure is unbounded in the x-y plane; 

hence a solution to the potentials is given by (see, for example, ref. 8, p. 145) 
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(2 T ,) 2 


dk^ dky 




(2tt)^ _ 


dk x dky 
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where 


^(kx’V 2 ) = A ( k x> k y) e 
F (k x ,k y ,z) = B (k x ,kyje 


i k z nz 


i k z n Z 




( 4 ) 


The roots of k z ^ chosen in accordance with the radiation condition at infinity are 

k Z H = f O 2 - k x 2 - ky 2 ( k 0 2 >k x 2+ k y 2 ) 


k z n " ^\j k x 2 + k y 2 " k Q 2 ( k 0 2 < k x 2 + k y 2 ) 




(5) 


If equations (3) are substituted into equations (1) and the boundary conditions are 
applied at z = z Q (i.e., continuity of the transforms of H x , Hy, E x , and Ey), the 
fI ( z o) G*( z o) 

ratios , and their derivatives can then be determined. By defining 

F n (z 0 ) G n (zo) 


f(z) 


_ F ( z ) and g(z) = ^ ( z - )- ; the resultant boundary conditions 


fII ( z o) 


G :I ( z o) 


are 


F ( z o) 


f'( Z o) = 


- d l " Fl (z) 


S( z o) 


"Gn ( z o) - 


dz [F n (z 0 ) 

g ’( z o)=i 


Z=Z r 


d 

G X (z) 

dz 

G U (Zo)_ 


= ikz 11 


_ ik n e i( z °) 
~ ]z ~TT~ 




L ]z=z. 


( 6 ) 


where the roots of k z 1 ^ are given by equations (5). Equations (3) also reduce equa- 
tions (2) to 
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f = 0 


(7a) 


d 2 f 

dz 2 


n e j( z ) 2 2 

k “ 1 _ k z _ k “ 

|K 0 e n K x K y 


d 2 g _ l de i( z ) dg . 
dz 2 6j(z) dz dz 


u 2i£_ 2 _ k 2 

*0 e Q *x K y 


g = 0 


(7b) 


where f(z) = and g(z) = 


pII ( z o) 


G n ( z o) 


Since equations (6) and (7) are identical in form to the plane-wave expressions 
given by Swift and Evans in reference 6, the procedure given therein may be used to 
determine g(z), f(z), and their derivatives at any point in the slab as a function of /3, 
where k Q 2 j3 2 = l^ 2 + k y 2 . The fields in the inhomogeneous slab are therefore functions 

only of the mode constants G n (z 0 ) and F n (z 0 ) which can be solved by applying the 
boundary conditions at z = 0. 

If it is assumed that the total aperture field is the same as the incident TEqj 
waveguide mode, then the transform pair at the aperture is 


^o wy 

E v = — ^ cos — — 
x a b 


E x = 27rbV 0 


sin^-g— jcos 


: x a\ kyb 

2 J COS 2 


¥)[■* - Ml ! 


(Ey=0) 

(Sy = o) 
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From equations (1) and (8), the boundary conditions at z = 0 therefore result in 
the following relationships: ^ 

= f it k y G n (z 0 ) 1 

% = 0 = n F - ja , M()ei (o) gf(0) j 


JI 


E- 


' x = _j V 11 ^ + Vmo^i(O) 
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= 27rbV, 
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/] 

sinMH cos 
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k x a 


^ - (kyb) 2 




which allows F n (z 0 ) and G^Zq) to be expressed in terms of V 0 , as follows: 


fII ( z o) - , i ky Wo 


sin 
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(k x 2 + ky 2 ) [M ^ - (kyb) 2 
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ii f (°) 
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(10a) 
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gII ( z o) 


, v • /W a \ / k yb\ 

cos ^ - 

^ p -(kyb^ 


k x 

iW 


(10b) 


Therefore, the transforms of all the field components at the aperture plane are known to 
within a constant Vo, which is implicit in the transform of E x . By substituting equa- 
tions (10) into equations (1), the transform of Hy at z = 0 may be evaluated in terms 
of E x . The result is 


Hi 


= 1 . Y 0 E x k p 2 j/M 2 e i(°) k Qg(°) / k y\ 2 [ ~f’(0) 

" ( k x 2 + ky 2 )Vo/ € ° l- g ' (0) J _ V k 0/ W 


e r (0) 

k 0 g(0) 

e 0 

_g’(0) _ 


(ID 


The complex conjugate of P at the aperture in terms of the field quantities in 
region I is given by 

OO 

p *'i 11 Ex* H y I<b! d y = I ~2 JT % 1<lk x dk y (12) 

aperture ' 


where the dual integral representation is a statement of Parseval's theorem (ref. 8, 
p. 458). 

The complex conjugate of P at the aperture in terms of the fields inside the 
waveguide is given by 


P .-I Y kol 2 (l - r) P */2 
2 01 a 2 (l + r)J _ a/2 


r*b/2 

-b/2 


cos 2 f dx dy = i Y q1 


2 (1 - D b 
(1 + T) 2a 


(13) 


From continuity of power across the aperture, equations (11) and (12) give 


-12- = y. =g. -jb. 

Y 0 1 y in J in (i + r) 


2a 1 

_1_ _1_ 

b Y 01 

|V Q | 2 (2tt) 2 


If E * 


” *H y I dk x dky 


(14) 


The substitution of equation (11) and the explicit form of E x , as given by equations (8), 
into equation (14), leads to the following expression for the input admittance: 
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where a transformation to polar coordinates has been made such that k x = k Q 0 cos a 
and k y = k Q /3 sin a. The values of g(|3,0), g'(ftO), f(/3,0), and f’(/3,0) are provided by 
the numerical integration of equations (7), subject to the initial conditions (6). The 
admittance has been broken up into TM and TE components to denote that y is 

derived from the potential A z and that y TE is derived from the potential A Z/ /e . If 
surface waves are excited, it is convenient to study y^M and y^ E separately. (See 
the appendix.) 

If the losses within the slab are moderate, the integrands of equations (15) are well 
behaved functions and the numerical evaluation of the input admittance is straightforward. 
However, if the losses are small, it is possible that the integrands may become large 
(or infinite if there is no loss) at discrete values of /3. In the limit of zero loss, these 
points become poles of order one, and their residues become a measure of the power 
which is confined within the slab as surface waves. The discussion of this surface-wave 
component is deferred to the appendix. 

NUMERICAL PARAMETERS 

The differential equations (7) were numerically solved as a function of /3 by using 
the Runge-Kutta method. These numbers were stored on tape and were then used as 
input information for the integration of equations (15). Approximately 250 values of (3 
were selected in the interval 0 ^ ^ 6. From the analysis of a similar problem (ref. 9), 

it was believed that the termination of the integration at j3 = 6 would yield a susceptance 
error of about 1 percent or less. To test this supposition, a numerical check against the 
methods described in reference 3 was initiated, and the results are discussed later in 
this section. Including the time needed to solve the differential equations (7), each admit- 
tance point required approximately 45 minutes of machine time. 
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The plasma-covered antenna was of immediate interest; therefore, the dielectric 
constant of the inhomogeneous slab was assumed to be given by 



which is based upon a point application of the Lorentz model (ref. 7, pp. 325-327) of wave 
interaction with electrons in the plasma. With the dielectric constant defined by equa- 
tion (16), the admittance of an aperture covered with a homogeneous plasma was computed 
as a function of (a>p/u>)2 5 0.7 with u/u = 0.4 and z Q = 3.5 cm. On the basis of an 
' i ' -2oi 3 z 0 

empirical criterion of e v > 0.1, where Op is the attenuation constant, it was 
believed that the admittance would be independent of thickness over the computational 

range. The criterion was, indeed, satisfied because the slab computations were within 

1 9 

1- percent of the corresponding half- space calculations performed by R. C. Rudduck. 

This agreement was considered to be a satisfactory check of the slab program. 

In order to observe the influence of an inhomogeneous plasma slab on the input 
admittance of the antenna, the series of electron density profiles shown in figure 2 were 
considered. Such profiles are typically produced whenever a plasma flows tangentially 
over a flat plate, as illustrated in the inset of figure 2. The variation of the electron 
density near the ground plane is due to viscous boundary- layer effects. 

With the aperture dimensions fixed at 0.4 X 0.9 inch (inside dimensions of a stan- 
dard RG 52/U waveguide) and a fixed propagating frequency of 10.0 GHz, admittance cal- 
culations were performed as a function of ^ P ?Q j = ^-97 X 10 _j ^ ^ for g^^ 0 f f^g 
distributions shown in figure 2. In this way, it was possible to systematically observe 
the changes which can occur in the input admittance as a result of changes both in the 
peak electron density and the extent of the boundary layer. 

The thickness of the slab was again chosen to be 3.5 cm, and calculations were 
performed over a range of parameters which satisfied the empirical criterion. 


iThis criterion was deduced from the inspection of similar calculations (ref. 9) 
pertaining to long slots on plasma- covered ground planes. 

^The calculations were made in fulfillment of the grant reported in reference 10 
but were not included in the final published report. 
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RESULTS 


The calculated input admittance for v/w = 0.4 is presented as a function of peak 
plasma frequency and boundary-layer thickness in figure 3(a). The curves are labeled 
"profile 2a," "profile 2b," and "profile 2c" to indicate that the computations apply to each 
of the three electron density distributions given in figure 2. Figure 3(a) shows that a 
substantial reduction of the susceptance can occur as the boundary- layer thickness grows 
from zero to a small fraction of a free-space wavelength. At ( w p,o/ w j = for 
example, the normalized input susceptance resulting from the homogeneous profile 2a is 
-3.37 as opposed to -1.90 from the inhomogeneous profile 2b; yet the half-maximum of 
the electron density of profile 2b occurs only Xq from the ground plane. A further, 
but less pronounced, decrease of the susceptance is noted as the boundary- layer thick- 
ness increases to ^ Xq (profile 2c). 

( Figure 3(a) also shows that the conductance is relatively insensitive to the detailed 

i distribution of the plasma if the boundary- layer thickness is ~ ^ Xq which indicates 
[ that homogeneous slab approximation adequately defines the conductance if the preceding 
/ condition is satisfied. 


The voltage reflection coefficient for v/u = 0.4 is given as a function of peak 
plasma frequency in figure 3(b). The results show that a definite lowering of the reflec- 
tion coefficient occurs as the boundary- layer thickness increases, which is consistent 


with Taylor's experiments (ref. 4). It is of further interest to note that the slope 

dir I /» D 0 \ 2 

— 1 — n decreases in the region is I I ^10 as the boundary layer grows. The 

d (“p.o/") 2 

determination of electron density within a given percentage of error therefore requires 
a more accurate measurement of d|r| as the boundary -layer thickness increases. In 

d|r| 

l-V. « 1 -I w. -i 4- n ' ' H y\V«n nn lM-frinivi a4-i i n ■nrvniil'nn/'l f /o r\ /■! -i i nl rtnfnrvvi 


the limit as 


density. 


'(>» /<■>)' 


— 0, phase information is required to deduce the electron 


The collision frequency ratio v/u was reduced to 0.06, and admittance computa- 
tions were performed by using profile 2c of figure 2. Calculations were initiated at 
( w p = 2*0 where the 3.5-cm slab is thick enough to appear as an infinite half space 

in accordance with the criterion e ^‘ a 9 z ‘° > o.l. The results of these admittance com- 
putations are shown in figure 4(a). For reference, homogeneous-half- space calculations 
(ref. 10) for v/w = 0.04 are also given. 


The behavior of the susceptance is similar to the higher collision results of fig- 
ure 3(a), which indicates that the susceptance is relatively independent of the collision 
frequency. On the other hand, the conductance decreases with a decrease in collision 
frequency. 


12 



It is interesting to note that a finite standoff distance of the plasma prevents gj n 
from approaching zero when the collision frequency is low. This observation is impor- 
tant, since this small increase in g^ n is the prime reason for the reduction in the 
reflection coefficient shown in figure 4(b). It is also interesting to note that |r| is 
relatively insensitive to changes in electron density when (oj p,o/ w ^ > ft * s there- 

fore concluded that phase information is required for the determination of electron den- 
sity when the collision frequency is low and the plasma is overdense. 

Finally, it should be noted that the integrand of y^M was inspected in order to 
note any surface-wave excitation. No evidence of surface waves was observed in the 
region 1 < /3 < 6. It is possible that a surface-wave pole occurs beyond the range of the 
numerical integration (i.e., /3 > 6); however, it is believed that such a surface wave 
would be weakly excited. Hence, the observations noted previously are not influenced by 
the question of surface waves propagating in the boundary layer. 

CONCLUDING REMARKS 

The input admittance expression of a rectangular waveguide-fed aperture antenna 
covered with an inhomogeneous plasma slab has been formulated, and numerical results 
have shown that a finite boundary layer over the ground plane can significantly influence 
the input admittance of the antenna. In particular, a substantial reduction of the induc- 
tive susceptance occurs as the plasma becomes more overdense and the boundary- layer 
thickness increases from zero to a small fraction of a free-space wavelength. A cor- 
responding decrease in the reflection coefficient occurs for a collision frequency ratio 
v/u> of 0.4. 

When the collision frequency ratio is reduced to approximately 0.06, a decrease in 
the reflection coefficient is also noted. For this low collision frequency ratio, the 
decrease in the reflection coefficient is primarily due to the fact that the boundary layer 
prevents the conductance from approaching zero. It is further noted that, when the 
plasma is overdense, the reflection coefficient is relatively insensitive to changes in the 
peak plasma frequency, and phase information is thereby required for plasma 
diagnostics. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 12, 1967, 

129-01-03-03-23. 
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APPENDIX 


ADMITTANCE EXPRESSIONS FOR HOMOGENEOUS DIELECTRICS 


If the dielectric constant is independent of position, the Helmholtz equation (7) may 
be solved to give 


(Al) 


for N z > 0, where 


and for N 2 < 0, 


, x i k z Iz 

-jk z z 1 

f(z) = Ae 

+ Be l 

, x jk z I z 

; -jkjz j 

g(z) = Ce 

+ De Z J 

T 

— = ±\/n 2 - j3 2 

^0 

= ± i\ 

/f3 2 - N 2 

k i 


k z 

h 3 2 + N 2 

*0 V 

A, B, C, and D in equ 


(n 2 > |3 2 )l 
(n 2 < /3 2 ) 


(A2a) 


(A2b) 


from the boundary conditions (6) at z = z Q , and explicit forms of 


mm 


k Q f(/3,0) 


kQg(/3,0) 


g'O,0) 

may then be substituted into the admittance expressions (15) to give 


and 


.TM 4 i( k 0 a )( k 0 b ) pr* 
Ym /Y a ^0 ^0 


cos 


j3k Q b 


sin a 


-|2r 


0i / Y 0 J 0 ^0 |*2 . ( k()b/3 gin a £ 


sm 


k n a^ 


cos 


I 


( k oa/3 


cos a 


N 


*• 


1 _ i N 2 \/l - § 2 tan k z \/n 2 - 
0 ° V __ 




^N 2 - ( 3 2 


r 

tan k. 


0 z oVN 2 - /3‘ 


2 + j N 2 lT g 2 


- > cos 2 o! j 8 dj3 do? 


\/n 2 - p* 


(A3) 
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TE 


-4j( V)( k o b ) r°°p 
Y 0l/ Y 0 ^0 ^0 


cos 


/3k n b 


sin a 
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7 r 2 - ^k Q bj3 sin otj 


sin 

/k 0 a /3 \ 

\ 2 COS a J 

s 

k Q a^ \ 

— — cos aj 


n 2 


^N 2 - /3 2 


X 


tan k Q z 0 ^N 2 - /3 2 + j -f=== 

' 2 _ o 2 


Vn 4 - /rj 


1 _ i ^ ~ ft - tan k n z 
0 


q \/n 2 - £ 2 


> sin 2 Q! j3 d/3 da (A4) 


The input admittance is the sum of the foregoing TE and TM components - 
that is, 

TE TM 

y in = y + y 


(A 5) 


If the slab is lossy, 




and 


]/n 2 - (3 2 =■ 4 


i 

€ I? r 


e 0 




’* fe ) 


1/2 




/3 2 -^W(^-/3 2 ) 2 + fi^) 


e o ■ yv e o 


1/2 


and no difficulty is pres ented i n the evaluation of equations (A3) and (A4), other than a 
proper root change of \ 1 - /3 2 , in accordance with equations (5). 

If the dielectric slab is nonlossy, the integration over /3 in equations (A3) and 
(A4) should be broken up in accordance with the root conditions in k^kQ and k^^ko, 
as summarized in tables I and II for the following cases: 

Case 1: N 2 > 1 
Case 2: 0 < N 2 < 1 
Case 3: N 2 < 0 
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rpP TM 

It is noted in the tables that poles occur in the integrands of y and y for 
N 2 > 1 and in the integrand of y™ for N 2 < -1. These poles are of order one and 
occur at those values of 0 = 0 n which are the roots of the following transcendental 
equations: 

TM component: 

tan k n z 0 ^N 2 - fi 2 - n2 ^ 2 ' 1 = 0 

TE component: 

1 + p 2 ' 1 tan k„z 0 ^Tg2 „ 0 

TM component: 

tanh k Q z 0 \/|N 2 | + 0 2 = ~ 1 

y[N 2 | + 0 2 

Since the integrand is infinite at 0 = 0 n , it follows that the numerical integration 
of yTE an d yTM should be terminated at 0 n - 6 and commenced again at 0 n + 6, 
where 5 is an arbitrary small number. This procedure also necessitates adding a 
residue contribution at 0 = 0 n . The poles are of order one; hence, the residue of y™ 
for N 2 > 1, for example, is given by 


(n 2 > l) 


V (A6a) 


(n 2 > l) 


J 


(n 2 < -l) 


(A6b) 


N 2 /3n| 


Res(y™) N2>1 = 


t2,//3 2 
1 

Jn 2 - 0 n 2 


N Fn - 1 ,n 2 

1 + —!■?. tan k 0 z o VN 2 - 0 n 2 




tan k 0 z o \p - 0 2 - N ^ 2 - 1 


i/n 2 - 0 2 _ 




P = Pn 


,cos a da 


(A7) 


where 1 < 0n < N, and A(o,0) for 0 = 0 n is defined in table I. 
From use of equations (A6a), it immediately follows that 


1 + N ~ 1 tan k 0 z o JN 2 - 0 n 2 = sec 2 k 0 z o ^N 2 - 0 R 2 

\|N 2 - 0 n 2 


(A8) 
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which simplifies the numerator of equation (A7). The derivative of the denominator of 
equation (A7) is given by 


dj3 


r 


'\|n 2 - /3 2 

t 

tan kQZ 0 \ji 


\/n 2 - (3 2 

0 


0=0n 


n 


f 2 -(Sn 2 _ 


tan k 0 z o \|N 2 - ;5 n 2 - 


N 


w 


^N 2 - 0 n 2 


0n 


kqZo sec 2 k 0 z o \(N 2 - /3 n 2 + N , = + 


2 N 2 ^ n 2 - 1 


$7^1 ' N2 ^n 2 


(A9) 


From equations (A6a) the first term on the right-hand side of equation (A9) is, by defini- 
tion, zero. Further use of equations (A6a) gives 


N 2 N^n 2 " 1 _ ( n2 ' k 0 z o \I n2 - 0n 2 

H" n o ' " 

(0n 2 " # 


n2 - 

It is now easy to show that 


P 2 ~ 0n 2 


(A10) 


d/3 


r 

r 

\|n 2 - j3 2 

tan k 0 z o \jl 

l 



= -0, 


Jn 2 - /3 2 


0=0n 


= -0 n k O z o sec\z 0 \lN 2 - |3 2 


1 + (n 2 - l) sin ( 2k o z o^N 2 - 0 n 2 ) 
K 2 - l) (2k 0 z o \/N 2 - /3 n 2 ) _ 


(All) 


and that the residue is 


Res 




2 >1 = 


N 2 ^ A^Q!,j3 n jcos 2 Q!dQ! 


k 0 z ° 

L 

(n 2 - 1) 

sin(2k 0 z o ^N 2 - 0 n 2 ]] 

i 

^n 2 - i 

(2k 0 z o \/N 2 - 0 n 2 ^ 



(A 12) 
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Therefore, the pole contribution is 



The right-hand side of equation (A13) is a real number, contributing only to the 
conductance, and is a measure of the amount of power which remains trapped in the 
dielectric sheet as a surface wave. 

The onset of a TM surface wave occurs when j3 = 1 and tan ]£qZ 0 ^N 2 - /3^ is a 
positive number, that is, where 

mvr = k Q z 0 V& 2 - 1 (m = 0, 1, 2, . . .) (A14) 

The total surface-wave contribution is given by 

TM V TM ,. ic , 

g s =Z g s,n < A15 ) 

n 

where the sum includes all the poles as determined from equation (A6). The total TM 
admittance is therefore given by 

y™=j^ dacos 2 ®^ A(c^/3)Fj(/3)j3 d/3 

+ P j^ N A(a,(3)F 2 (0)f3 d/3 + J Afa/3)F 3 (/3)/3 d/s! g™ (A16) 

N J n 

where the explicit functional forms of Fj, F 2 , and F 3 are given in table I. 

The surface-wave conductance for y^' E |j^2>i anc * ^^In 2 <0 iS sim ^ ar ^y 
derived, and explicit forms are given in table in. 

If a large number of surface-wave poles exist in the region 1 < /3 < N, F 2 (jS) 
becomes a rapidly varying function. As a consequence, it may be difficult to evaluate the 
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principal part of the integral. If such difficulties are encountered, it would be wise to 
convert the integration over the real axis to one along the branch cut, as indicated in 
figure 5.^ 


If the medium is inhomogeneous and nonlossy, one should inspect g'(/3,0) and 
f(/3,0). If they go to zero, a surface wave has been excited, and one should numerically 


compute the values of 
residue contribution. 


k Qg(/3n>°) 



and 



and numerically enter the 


^ It should be noted that only the forward wave poles are included when the contour 
is closed in the upper half plane. When the contour is closed in the lower half plane, the 
backward wave poles must be included. These backward wave poles and the branch line 
which appear in the lower half plane are not shown in figure 5. 
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TABLE I.- TM COMPONENT OF ADMITTANCE EXPRESSION 


Preferred form for the numerical 
evaluation of y™ 


F lO) 




Case 1: y™ = f dor cos 2 <*[f A(a, p)Fi(0)p dft 

N 2 > 1 0 l: 0 < 

r N N 

+ J i A(a,p)F 2 mW _ . 


Comments 


l^j^^tank 0 z o lfi^7 


r ‘ 2 (~2 , ~ \ " 9 i “) (1) Poles of F 2 (p) occur where 

{2} principal parts of the inte- 


ftan h 


Vo^ 2 - ^ - ^ - |2 - r r l $ - N2 tanh N 2 + 

_ ^2 l u n; o 


Case 2: y™=f da cos 2 a f' N Afo^F^ d/3 

o < n 2 < l 0 irO < 

pi 

+ J M A fo0F 2 <00d0 

+ £ AteflFjOMd/jj 


grals which contain F 2 (0) 
and Fg (/3) contribute only 
to the susceptance. 


1 * 3 tan kqZ Jn 2 - ft 2 j ~N 2 j 

VN 2 - ff 2 J 


tanh k 0 z o ] h 2 ~ N 2 -N 2 [ 1 + > 


tan k Q z 0 VN* - /3 2 + j 


r 2 . *2 + 3 - N 2 Vl - £ 2 


JZ3 


in*- a 2 J 


r tanh kqzjff - N' 


tanh kgZ ( 


- N 2 - j 


Case 3: Y™ * f da eos 2 a[\ d0 |n 2 | 1 + j ^ tanh k Q z ^[n^T? 

N 2 <0 0 irO |_ ^)n 2 j + j9 2 ° 

^ AM)F3 H ^L Vo ii^ +i iNilOl 

u ,/lT5T 


I N 2 1 + ^ 2 


A(a,|3) 


_ 4 i(V)M 


y 0l/ Y 0 j^ 2 - (^(3 Sin a) 


nw , i?r /v» \? 

\~2~ sin aj sinl -Tp cos a j 


(-T- cos “) _ 


^ 2 -n 2 


- N 2 tanh k 0 z o \j & 2 - N 2 + N 2 ^ - l j 
u »nr o 


V - N 2 


(1) No poles occur anywhere in 
the range of integration. 
L. (2) The integral which contains 
F3C/3) contributes only to 
the susceptance. 


(n 2 | 


]/|n 2 | + 0 2 ! 


(1) Pole of Fg(jS) occurs where 
1 - tanh vi^U? 2 ] tanh Vol/i^TT?- i*d!fiU 

(U 2 | > 1) 


tanh kg! 


)Zo nPu 7-!^S 


i!In 2 | + /3 2 


(2) The integral which contains 
F^W) contributes only to 
the susceptance. 


Note that 


sin(^ cos a) 


(V? y 

{—cosaj 


a 1 for j9 cos a ~ 0 and 




n 2 - (k Q b0 sin aj 


g for k 0 b @ sin a ~ 5r - 

16ir 





TABLE IL- TE COMPONENT OF ADMITTANCE EXPRESSION 


Preferred form for the numerical 
evaluation of y^ E 


F 1 (/3) 


f 2 (0) 


f 3 M 


j n? — n 

-\jp 2 - N 2 

tanh k n Zn\fs 2 - N 2 + " 1 

[ fa - M*l 

i J 

r 

^i-taihk^o^-N 2 


Case 1: y 1 

N 2 > 1 


X dasl " 2 “to 

r N 

+ A(a J p)F 2 (p)p dp 

J 1 

1 

+ J A(a,j3)F3(0)0d/9 


dp 


i^ 2 - p 2 


|tan , 0 z> 2_ - p 2 + j - jlK" 


fe 2 - g 2 l 


\[n^ 


1 _ j ^ P tan kQZo^N 2 ”- j3 2 

if 2 - p 2 


1 + JZ n 

'In 2 - p 2 


1 Vo'In 2 - P 2 - f= ~ 2 1 

'/n 2 - J3 2 j 


1 k o z o' 


.Vn 2 - 0 2 


(1) Poles of F 2 (/3) occur where 


(2) Principal parts of the inte- 
grals which contain F 2 (/3) 
and F 3 (/3) contribute only 
to the susceptance. 


Case 2: ,y T 

0 < N 2 < 1 


9 fpN 

= -\ dosing A(o<0)F X (P)P dp 

J ° - J ° 

+ y A{a,p)F 2 mAI> 

”» 1 1-) ^-l^- tank 0 z o '|N 2 - g 2 

+ j i A(a,/3)F 3 (0)0d0 \|n 2 - p2 


: Q z 0 ^N 2 - p 2 + j - N 2 -tarih k^e 2 - N 2 + j ■j 1 r& 

u .Co n u .To 0 


fa 2 - a 2 


1 - j j 1 -" ^ tanh y-z 
(p 2 - N 2 


'Jp 2 - N 2 





tanh VqZ 0 \jp 2 - N 2 + ^ ~ * 

ip 2 - N 2 


P 2 ~ 1 tanh k 0 z o \|/3 2 - N 2 


(1) No poles occur anywhere in 

the range of integration, 

(2) The integral which contains 

F 3 (/3) contributes only to 
the susceptance. 


Case 3: y TE = -£ da sin 2 a J A(a,0)F ^p)p dp ^In 2 | + p* -tanh k Q z 0 \/lN 2 | + p 2 + ] ^ ~ ? 

N 2 <0 0 L0 [ \11n 2 1 + 


•■j* A(a,p)F 3 (p)p 


dp 


1 - j ^ tanh k 0 z o \fjN 2 |T^ 2 


/In 2 I + p 2 


tanh 1 CqZ 0 \ 




v/In 2 I +• p 2 _ 


jfif. - . i , 


/|n 2 | + p 2 


tanh k n z n v|N 1 + f. 


(1) No poles occur anywhere in 

the range of integration. 

(2) The integral which contains 

F 3 (/3) contributes only to 
the susceptance. 



TAkLE m.- TE AND TM SURFACE- WAVE CONDUCTANCE EXPRESSIONS 


Expression for surface-wave conductance 


Transcendental equation for the 
location of the roots /3 n 


Case 1: 



(1) The sum on n of g s includes all finite 

roots of the transcendental equation. 

(2) The onset of a surface wave occurs at thick- 

nesses in wavelengths ~ * 0 , 1 . 

. , * 

(IN 2 - 1 2\/n 2 - l’ 


(1) The sum on n of g s includes all finite 



to 


CO 
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Profile 2 c 


/ 



v/co = 0.04 for homogeneous half space 
v/u> = 0.06 for profile 2 c 


I I J.J.1.1 
10 ° 

%„« 2 


J L 


10 J 


(b) Voltage reflection coefficient. 
Figure 4.- Concluded. 
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